Spin-dependent Hedin's equations 
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Hedin's equations for the electron self-energy and the vertex were originally derived for a many- 
electron system with Coulomb interaction. In recent years it has been increasingly recognized 
that spin interactions can play a major role in determining physical properties of systems such as 
nanoscale magnets or of interfaces and surfaces. We derive a generalized set of Hedin's equations for 
quantum many-body systems containing spin interactions, e.g. spin-orbit and spin-spin interactions. 
The corresponding spin-dependent GW approximation is constructed. 

PACS numbers: 



An increasing number of modern applications of mate- 
rials science imply spin degrees of freedom. Spintronics 
devices, colossal magneto resistance materials or mag- 
netic impurities in semiconductors are important exam- 
ples, not only due to their technological interest but also 
owing to the fundamental questions they raise. In fact, 
a detailed understanding of such systems requires a the- 
oretical description of the interplay of spin, charge and 
orbital degrees of freedom, and relies crucially on a re- 
liable description of electronic interactions. In striking 
contrast to these scientific and technological needs, the 
actual progress in developing techniques for describing 
interaction effects beyond density functional theory [l[ 
in an explicitly spin-dependent manner has been slow. 
In particular, one of the most important first principles 
many-body methods used nowadays for describing elec- 
tronic excitations, the so-called GW method, still awaits 
its extension to spin-dependent interactions. 

In 1965, Hedin derived a closed set of equations for the 
electronic Green's function and self-energy, the screened 
Coulomb interaction and the polarization of a solid 0. 
The equations provide an iterative scheme for an expan- 
sion of the self-energy in powers of the screened inter- 
action. Although there is no guarantee that the iter- 
ative procedure is absolutely convergent, as originally 
envisioned by Hedin, it nevertheless provides a rigorous 
basis for studying the self-energy of real materials from 
first-principles. In particular, the lowest order approxi- 
mation leads to the GW approximation (GWA) 01 
which has proven very successful in studying one-particle 
excitation energies of real materials entirely from first- 
principles. 

The original Hedin's equations were derived for a 
many-body Hamiltonian with a Coulomb interaction 
only, without the possibility of having an explicit spin- 
dependent interaction. The GWA derived from them 
has met tremendous progress in describing the electronic 
structure and excitation spectra of semiconductor-based 
systems, transition metals and oxides. Due to their re- 
striction to spin-independent interactions it has however 
not been useful for studying correlation effects arising 



from spin-dependent interactions in the classes of ma- 
terials mentioned above. These interactions are crucial, 
despite the tiny energy scales associated with them. To 
cite an example, the conduction band spin splitting in 
zinc blende semiconductors arising from spin-orbit cou- 
pling is only tens of meV and yet it is important for 
applications in spintronics since it determines the spin 
lifetimes as well as inducing spin current in the absence 
of a magnetic field, the so-called Rashba effect [H . 

The aim of the present work is to present general- 
ized Hedin's equations, suitable for dealing with explic- 
itly spin-dependent interactions. These interactions may 
arise from relativistic effects, such as spin-orbit coupling, 
or from an external perturbation like in the case of a 
magnetic impurity in a semiconductor. With a spin- 
dependent interaction it is possible to relate the two- 
particle Green's function to functional derivatives of the 
one-particle Green's function with respect to an electric 
and a magnetic field. The first order term in the screened 
interaction of the resulting spin-dependent Hedin's equa- 
tions yields a generalization of the GW approximation. 

The Hamiltonian with a spin-dependent two-particle 
interaction is given by 



d'r ^jj^ (r)/io(r)V„(r) 



(1) 



where ho is the one-particle Hamiltonian. The identity 
of the particles implies that the second term is invariant 
under particle interchange: ?;K'y/3,,(r, r') w^^K7(r',r) 
0. For this general Hamiltonian, we will derive the fol- 
lowing set of generalized Hedin's equations, relating the 
electronic self-energy E to the Green's function G and the 
screened interaction W, using the polarization P and the 
vertex function A: 



S„/3(l, 2) = -<Ti^a„(l, 4)A:J^(4, 2, 5)Wjii5, 1), (2) 
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The polarization then becomes 



Wij{\, 2) = vij(l, 2) + viK{l,i)PKL(i, 4.)Wlj{4., 2), 



(3) 



Pzj(l, 2) =ai^a/3^(l,3)A4(3, 4, 2)0,44,1+), (4) 



Ai.(l,2,3) = 5(l-2)<5(2-3)ai 



a/3 



+ ^^^^e,v(4, 6)A;,,,(6, 7, 3)g.,(7, 5). (5) 

Here, ct*, i=x,y,z, are the Pauh spin matrices and cr" 
is defined to be a 2x2 unit matrix. Capital letter indices 
run over 0,x,y,z, while Greek letters take the values ±1. 
We use the following common shorthand notation: (xr) 
is represented by a number, repeated indices are summed 
and repeated variables represented by numbers are inte- 
grated, unless they appear on both sides of the equation. 
For clarity, we further adopt a notation that quantities 
with subscript denoted by capital letter do not depend on 
spin. The spin-dependent interaction has been expanded 
in the Pauli and unit matrices as 



ya,^^(l,2) = ai vij{l,2)a 



(6) 



The spin-Hcdin equations can e.g. be applied to Hamil- 
tonians containing interactions of the following form. 



Wa^/3r,(r, r') = 



aV^^/|r-r'|, 
CT^^Jy(r,r')cr^,, 



(7) 



where the first is the usual Coulomb interaction, the sec- 
ond a spin-spin interaction and the third a spin-orbit 
interaction, which contains the angular momentum op- 
erator. 

The general structure of the set of equations ([1]) - © 
is remarkably close to the usual Hedin's equations. How- 
ever, the self-energy now depends on the spin variable, 
and the screened interaction W as well as the polarization 
function P acquire a matrix form expressing an interplay 
between the charge and spin channels. Thus, a polariza- 
tion function Pqi: for example, describes a charge density 
response of the system with respect to a perturbing mag- 
netic field in the i-direction. When P is used in ([3]) the 
charge channel of the screened interaction experiences 
the effects of the spin interactions and vice versa. 

Most importantly, one can construct from these equa- 
tions a spin-dependent generalization of Hedin's GWA, 
by approximating the vertex functions by 



Pjj{l, 2) = ai^gp-^il, 2)a4a,a(2, 1+) (9) 
yielding the self-energy 



(l,2) = -f7^,e,7(l,2)a^^l^j,(2,l) 



(10) 



A^4l,2,3) = ^(l-2)5(2-3)ai^ 



(8) 



Before proceeding further, let us interpret the meaning 
of the resulting polarization and the spin-dependent self- 
energy. Consider first the case I — J — Q (charge 
channel) giving Poo(l,2) = 5a^(l, 2)(?.^a(2, 1+) and 
E„/3(l,2) = -^a/3(l,2)M^oo(2,l). If the Green function 
is diagonal in spin space we recover the polarization and 
the self-energy in the original Hedin's equations. How- 
ever, for a system with an existing spin structure, such 
as a non-colinear spin, the Green function possesses non- 
diagonal spin components. This case is a generaliza- 
tion of the original Hedin's equations to spin-dependent 
Green's function and self-energy with purely Coulombic 
interaction. It emerges naturally in the present formula- 
tion as a special case where spin interactions are absent. 

Let us now consider the case when the interaction is 
spin dependent, which may arise from purely spin-spin 
interaction or spin-orbit coupling, among other possibili- 
ties. The exchange-correlation effects on the Green func- 
tion of up spin is illustrated in Fig. 1 (lower panel). A 
particle of up spin enters the self-energy Upon 
entering the self-energy the electron spin is flipped to 
down spin by a spin operator cr||and a magnon repre- 
sented by Wij is emitted. Upon leaving the self-energy 
the spin operator cr|| causes the electron to reabsorb the 
magnon and return to its original up spin configuration. 
This process is analogous to the original Hedin's GWA 
whereby an electron emits and absorbs a plasmon but 
without the possibility of spin flip (upper panel). 

The derivation of the generalized Hedin's equations 
closely follows Hedin's original work, using Schwinger's 
functional derivative technique. Since many problems re- 
lated to spin degrees of freedom involves temperature, we 
work in the finite-temperature formalism but of course 
the zero-temperature version readily follows. Using the 
Heisenberg equations of motion we obtain 



— + /io(x) - ^ 5q/3(xt, x't') 

OT J 

- j d^r v^^ocriiY, x)a^^^^„(xT, x'r', rr, vt+) 

= -<5,/3J(x - x')(5(t - r'). (11) 



To utilize the Schwinger functional derivative technique 
we work in the Dirac or interacting representation and 
define the Green functions as follows. 
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that 





FIG. 1: Diagrams for the self-energy within the spin- 
dependent GWA, compared with the conventional GWA see 
text 



(16) 



gg^Jl,2,3,3+) <(17) 



This important relation allows us to replace the two- 
particle Green's function by a functional derivative of the 
one-particle Green's function with respect to the applied 
electric and magnetic field. Using the above relation we 
define the mass operator A4 as follows. 



a7] 



The charge and spin density is given by 
Using the identity 



,(18) 



(19) 



r [5V'Ji)V'^(2)] 



S 



(12) 



,2, (r [5V'J1)V,(3)V'^(4)^A^(2)] 

eiX,(l,2,3,4) -TP 

'S 



(13) 



where T is the imaginary-time-ordering operator and 

r/9 



S = T cxp 



dr 



(14) 



In the original Hedin's derivation a probing electric 
field was applied to obtain an equation relating the two- 
particle Green's function to the functional derivative of 
the one-particle Green's function with respect to the ap- 
plied electric field. Since we have spin-dependent in- 
teractions we have to consider not only a probing elec- 
tric field, ipQ{rT), but also a probing magnetic field, 
(Pj(rr), i = x,y,z. The electric and magnetic fields are 
given by 



)= d^r ifi{vT)^^{vT)aipi;p{YT). (15) 



Noting that 5S /5ipi{l) — —T[S(7 (1)] it can be shown 



{5g/5^) = ~g{5g-^/5ip)g, 



(20) 



multiplying by g^^ on the right and taking out the 
Hartree-like term (the first term) we obtain from (HI 
an exact expression for the self-energy: 



i:^i3{i,2) = vij{i,i)cj'g^^{iA)- 



^5,-/3' (4, 2) 



As pointed out in 0] in the absence of an explicit two- 
particle spin interaction, only / = J = components sur- 
vive and the self-energy depends only on variation with 
respect to the electric field (pg. Application of a magnetic 
field introduces a one-particle term in the Hamiltonian 
but it docs not affect the self-energy. 

Next, we identify the generalized Hartree potential 



Vj"{l)^vi.j{l,i)pj{?,) 



(22) 



and the total field (f>/ = (^j + V/^. The vertices can then 
be defined as follows: 



AiMl'2,3) 



^$,(3) ■ 



(23) 



As usual, the dielectric function is related to the deriva- 
tive of the total field with respect to the applied one 



<5$j(l) 
5vj{2) 



{l,2)^5{l-2)5ij + 



S^j{2) ' 



(24) 
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Thus, using the chain rule 



A;^^(l,2,4)e7;(4,3), (25) 



and inserting this equation into (PT|) we obtain Eq. ^ 
with Vl^/j(l,2) = e7f^(l,3)i;/fj(3,2). The charge and 
spin vertex equations are given by 



Ai,(l, 2, 3) = 5{l 2)5(2 - 3)ai, + ^||^ 



^Sa/3(1,2) 

(5^7^,(4,5) 



g^^-(4,6)A^,j6,7,3)g,^(7,5)i 

(26) 



where we have made use of the chain rule, the identity 
and the definition of the vertex in 
We now establish the equation for the screened inter- 
action. The polarization is given by the variation of the 
density with respect to the total field. Using ([TO]) . (PO)) . 
and 



Green's function, when an iterative scheme is used. For 
example, if we consider an Ising interaction and choose a 
starting Green's function corresponding to the mean-field 
solution, solving the magnetic Hedin's equations itera- 
tivcly will not generate other spin configurations because 
there is no spin-fiip term in the Ising model. More explic- 
itly, the polarization in ([9]) with the mean-field Green's 
function will only possess the z-component Pzz , which in 
turns will only generate Wzz- Thus, the self-energy in 
pop will not have off-diagonal components in spin space 
that would allow for spin-flip processes yielding spin con- 
figurations other than the mean-field one. To solve this 
problem, one must start with a starting Green's function 
that corresponds to a multi configuration. This prob- 
lem does not arise in the general spin spiral case, where 
the spin at each point may have different magnitude and 
direction. The non-interacting Green's function Q 

Qap (r, I- ; = 2^ — -—f — (30) 

n 

possesses non-diagonal components in the spin space 
which contribute to the charge channel of the polariza- 
tion. 



5*7(2) 

, 50^^(1,1+) 
'^"^ <5$,7(2) 

= -C^a/3^/37(l'3)-^^-^^^a(4, 1 ) 

= ai^a/37(l, 3)A4(3, 4, 2)g,„(4, 1+). (27) 
Expressing the dielectric function (|24p as 

<5Ff (1) 6puCi) ^*k(4) 



£7;(1, 2) = 5(1-2)57,7 + 



5p,,7(3) 5$7f (4) 5^,7(2) 



= 5(1 - 2)57,7 + i;7A7(l,3)PAf7f(3,4)£^i^(4,2) 

(28) 

and multiplying by the bare interaction, we obtain Eq. 
([3]). The response function, R = 5p/5(p, satisfies a similar 
equation: 

i?7,7(l, 2) = P7,/(l, 2) + P/Tf (1, 3)f7fL(3, 4)i?7f j(4, 2) 

(29) 

which follows from 5p/5(f = 5/9/5$ x 5<i>/5(y3 = Pe ^. 
The above three equations for e^^, W , and R are equiva- 
lent. This concludes our derivation of the spin-dependent 
Hedin's equations. 

In the following we first comment on some practical as- 
pects of the generalized GWA; we then conclude by dis- 
cussing possible applications and further developments. 

We note that while the magnetic Hedin's equations 
are exact, the solution depends crucially on the starting 



Poo(l,2) = alpg0^{l,2)a%^g^^{2A+) (31) 
= ^,^(l,2)g^„(2,l+) 

and the self-energy naturally acquires non-diagonal com- 
ponents in the spin space. 



Sf7(l,2) = -e„^(l,2)T4^oo(2,l). 



(32) 



Using the self-energy one solves the quasiparticle equa- 
tion 



TT 



S|'|^(-E'n 



(33) 

which contains the effects of the non-diagonal spin com- 
ponents of the self-energy. 

A non-interacting Green's function constructed from a 
Hamiltonian containing a spin-orbit interaction will also 
have non-diagonal spin components. A similar formu- 
lation as above can be applied to this situation. The 
screened interaction will contain charge-spin coupling 
components and the self-energy will take a more general 
form as given in PU)) . Correlation effects on orbital mo- 
ments and spin densities can then be accessed including 
life-time effects. 

In conclusion, we have derived the spin-dependent gen- 
eralization of the original set of Hedin's equations. Their 
first order term (in the spin-dependent screened inter- 
action) leads to a spin-dependent GW approximation. 
These equations allow for a truly first-principles study 
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of a wide range of problems where correlation effects in- 
duced by spin interactions play a crucial role in determin- 
ing physical properties. Applications to nanoscale mag- 
netic systems ranging from quantum dots, quantum wires 
to impurities or nanoparticles, as well as to filrns, surface 
and interface problems now come into reach [i,[lo|. Also, 
developments at the interface of many-body perturbation 
theory and time-dependent density functional theory as 
proposed in ll|, 112] are readily generalized e.g. to rela- 
tivistic interactions. The same is true for combined GW 
and dynamical mean field techniques [l3| . 
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